General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



118 5- 11425 


(£85** 10012 MASA-C1- 166569) MEflOTE SEISING 
IBVlLBD to fOBEST HESCOBCBS (Instituto de 
Pesquisas Espaciais, Sao Jose) 27 p 
NC AO 3/ At AQ1 Onclas 

G3/43 00012 



SECRETARIA DE PLANEJAMENTO DA PRESIDENCY DA REPUBUCA 
CONSELMO NACIONAL DE DESENV0LVIME4T0 CIENTIFICO E TECNOLOGlCO 




F85-10012 


RECEIVED BY 
NASA STI FACILTTY 
DATS 

I ill 111 ' ■■ 

DCAF Na 




processed by 

[t^NASA STI FACILITY 
□ ESA-SDS □ AlAA 


W INSTITUTO DE PISPUIGAS ESPACIAIS 


1. Publicagao n9 


2 . Versao 

3. Data 

5. Distribuigao 

INPE-3 1 6 7-PRE/53 7 



Julho , 1984 

□ Interna 0 Externa 

4. Or i gem 

DSR 

Programa 

RECSAT 


□ Restrit.a 


6. Palavras chaves - selecionadas pelo(s) autor(es) 

SENSORES REMOTOS 

REFLORESTAMENTO 

MAPEAMENTO 


7. C.D.U.: 528. 711.7:634.0.5 


8. Titulo 


INPE-3 1 6 7-PRE/53 7 


SENSORES REMOTOS EM LEVANTAMENTOS FLORESTAIS 


9. Autoria pecb } o Hernandez Filho 


Assinatura responsavel 








10. Paginas: 26 


11. Oltima pagina: 22 


12. Revisada por 


Antonio Tebaldi Tardin 


13. Autorizada por 


e Jesus Parade 
Dire tor Geral 


14. Resumo/Notas 

0 objetivo deste trabalho e apresentar o desenvolvimento de me_ 
todologias para classificar areas reflorestadas utilizando dados de sensoria 
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gao da idade (de: 0 a 1 ano } 1 a 2 anos t 2 a 3 anos s 3 a 4 anos t 4 a 5 anos 
e 5 a 6 anos) j da heterogeneidade do plantio (de: 0 a 20%, 20 a 40%, 40 a 
60%, 60 a 80% e 80 a 100%). ds diferengas relativas entre as areas com flo_ 
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dos com a finalidade de estimar o volume de madeira atraves de amostragem,em 
multiplo estagio , com probabilidade vaidada. 
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REMOTE SENSING APPLIED TO FOREST RESOURCES 


Pedro Hernandez Filho 


RESUMO 


0 objetivo deste trabalho e apresentar odesenvol 

vimento de metodologias para classificar areas ref lorestadas 

utilizando dados de sensoriamento remoto. Urn estudo preliminar 

foi conduzido no nordeste do Estado de Sao Paulo em 1978.0s re 

f lorestamentos de Pinus spp e Eucalyptus spp foram mapeados con 

base nas caracteristicas espectrais, espaciais e temporals das 

imagens do LANDSAT. Posteriormente , um estudo mais detalhado 

foi realizado no Estado do Mato Grosso do Sul. As areas reflo 

restadas foram mapeadas em fungao da idade (de: 0 a 1 ano, 1 a 

% 

2 anos, 2 a 3 anos, 3 a 4 anos, 4 a 5 anos e 5 a 6 anos) e da 
heterogeneidade do plantio (de: 0 a 20%, 20 a 40%, 40 a 60%, 60 
a 80% e 80 a 100%). As diferengas relativas entre as areas com 
floresta artificial, estimadas a partir dos dados do satelite 
LANDSAT e informagoes de campo, variaram de - 8,72 a + 9,49%. 
Atualmente, desenvolvem-se estudos com a f inalidade de estimar 
o volume de madeira atraves de amostragem, em nultiplo estagio, 
com probabilidade variada. 


ABSTRACT 


The objective of this study is to present the 
development of methodologies to classify reforested areas 
using remotely sensed data. A preliminary study was carried 
out in northeast of the Sao Paulo State in 1978. The reforested 
areas of Pinus spp and Eucalyptus spp were based on the 
spectral, spatial and temporal characteristics of the LANDSAT 
imagery. Afterwards, a more detailed study was carried out in 
the Mato Grosso do Sul State. The reforested areas were mapped 
in functions of the age (from: 0 to 1 year, 1 to 2 years, 2 to 
3 years, 3 to 4 years, 4 to 5 years and 5 to 6 years) and of 
the heterogeneity stand (from: 0 to 20%, 20 to 40%, 40 to 60%, 
60 to 30% and 80 to 100%). The relatives differences between 
the artificial forest areas, estimated from LANDSAT data and 
ground information, varied from -8,72 to +9,49%. Presently, 
studies are carried out with the aim to estimate forest volume 
through a multistage sampling technique, with probability 
proportional to size. 
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1 - INTROPUCAO 


0 Instituto de Pesquisas Espaciais (INPE) tem de 
senvolvido estudos aplicados a area florestal atraves da utiljL 
zagao da tecnica de sensoriamento remoto desde o ano de 1977, 
em convenio com o Instituto Brasileiro 'de Desenvolvimento Flo 
restal (IBDF) , com enfase em imagens dos satelites da serie 
LANDSAT. 


Com o intuito de desenvolver metodologias de sen 
soriamento remoto para avaliagao de areas ref lorestadas , em 1978 
foi escolhida uma area piloto, situada a nordeste do Estado de 
Sao Paulo, que abrange os municipios de Mogi-Guagu, Aguai,Casa 
Branca, Tambau, Santa Rita do Passa Quatro, Santa Rosa do Viter 
'jo, Sao Simao, Luis Antonio, Cravinhas, Serra Azul, Ribeirao 
Preto e Altinopolis (Hernandez Filho et alii, 1978). Posterior 
mente, o estudo foi estendido ao Estado do Mato Grosso do Sul, 
cuja escolha e justificada pela presenga de grandes concentra 
goes de areas ref lorestadas . Foram mapeados ref lorestamentos ao 
nivel de genero, idade e grau de heterogeneidade . Os resultados 
destes estudos estao documentados nos relatorios de Hernandez 


Filho, 1980, Shimabukuro et alii, 1981, e Hernandez Filho et 
alii, 1982. 


Atualmente, esta se avaliando o volume demadeira 
das areas ref lorestadas com base no inventario florestal emmul 
tiplo estagio, conforme a metodologia proposta por Langley 
(1975). Basicamente, para a aplicagao do metodo de inventario 
florestal, utilizam-se simultaneamente os dados contidos em ima 
gens orbitais e fotc^rafias aereas, complementadas com informa 
goes He campo. Este metodo de inventario florestal foi delinea 
do para gerar informagoes dos recursos florestais ao nlvel re 
gional, visando principalmente levantar dados para urn plane j a 
mento nacional de locagao de recursos e fornecer subsldios pa 
ra orientar o estabeleci mento de uma polltica florestal de urn 
pals. 

2 - DADOS DO SA T fiLTTE LANDSAT 

0 sistema LANDSAT, de visao sinoptica e repetiti 
va, foi desenvolvido para permitir a aquisigao de dados, comba 
se nas caracterJ sticas espectrais e temporais , sobre a super fl 
cie da Terra. 


Os tres primeiros satelites da serie LANDSAT, lan 
gados em 1972, 1975 e 1978, conduzem basicamente dois tipos de 
sensores: urn imageador multiespectral de varredura (MSS) de 4 
canais (5 no LANDSAT 3) e urn sistema de televisao (RBV) de 3 ca 
nais (1 no LANDSAT 3) . 




Os satelites tern orbitas quase-polar, circular, 
geocentrico e heliosslncrona, com perlodos de 103 minutos. As 
sim sendo, diariamente sao percorridas 14 orbitas, o que permi 
te um ciclo repetitive a cada 18 dias, ou 251 revolugoes, para 
que o satelite volte a imagear uma mesma superflcie da Terra. 

A Figura 1 mostra o esquema de passagens diarias 

do LANDSAT. 



Fig. 1 - Esquema geral das passagens diarias do LANDSAT. 


A configuragao do LANDSAT e mostrada na Figura 2. 
A imagem MSS corresponde a uma area aproximada de 35.000 km 2 , 
de 185 km por 185 km. 








Fig. 2 - Configuragao geral do LANDSAT. 


A radiagao coletada pelo satelite da superficie 
da Terra e decomposta, por meio de filtros, era 4 bandas espec 
trais ou canals, conforme a Tabela 1. 




TABELA 1 


FAIXAS ESPECTRAIS DO MSS DO LANDSAT 1, 2, 3 


CANAL 

MSS-BANDA ESPECTRAL 

4 

de 0,5 a 0,6 pm (Verde-laranja) 

5 

de 0,5 a 0,7 pm (Laranja-vermelho) 

6 

de 0,7 a 0,8 pm (Vermelho-infravermelho proximo) 

7 

de 0,8 a 1,1 pm (Infr^verm^lho proximo) 

8 

de 10,4a 12,6pm (L3) (Infravermelho termal) 


0 canal 8 e senslvel a temperaturas de -13°C a 
67°C com uma resolugao de 240 pm e senslbllldade de 1,5°C e so 
mente fpi concebido no LANDSAT 3. 

A forma de imageamento e mostrada na Figura 3. 

Os satelites 4 e 5 foram langados em 1982 e 1984 
respectivamente, conduzindo basicamente urn imageador multiespec 
tral (MSS) e urn mapeador tematico (TM) . 0 MSS apresenta as mes 
mas falxas espectrais que as de satelites anteriores. Asbandas 
espectrais do TM estao apresentadas na Tabela 2. 
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Fig. 3 - Forma de imageamento do sensor MSS 



TABELA 2 


BANDAS ESFECTRAIS DO TM DOS LANDSAT 4 E 5 


CANAL 

TM-BANDA ESPECTRAL 

1 

0,45 

a 

0,52 pm 

2 

0,52 

a 

0,60 ym 

3 

0,63 

a 

0,69 ym 

4 

0,76 

a 

0,90 ym 

5 

1,55 

a 

1,75 ym 

6 

o 

o 

a 

12,50 ym 

7 

• 2,08 

a 

2,35 ym 


3 - Area de estudo 

Para o desenvolvimento do trabalho foi escolhido 
o Estado do Mato Grosso do Sul, situado na regiao Centro-Oeste 
do Brasil (Figura 4). A concentragao dos reflorestamentos esta 
localizada no Distrito Florestal, em torno da rodovia R-262 e 
da Ferrovia Noroeste do Brasil, entre as cidades de Campo Gran 
de e Tres Lagoas. 
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Fig. 4 - Locdlizagao da area de estudo - Estado do Mato 
Grosso do Sul. 


As imagens do satelite LANDSAT que cobrem as areas 
de ocorrencia de ref lorestamento sao: 220/26 e 27, 234/26 e 27 
e 248/25, 26, 27 e 28, segundo o Sistema de Referenda do Bra 
sil para as orbitas/pontos dos tres primeiros satelites da se 
rie LANDSAT (Figura 5). 
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Fig. 5 - Imagens que cobrem a regiao de estudo. 


4 - LEVANTAMENTO DAS AREAS RFFLORESTADAS - PRIME IRA FASE 


No ano de 1979, foram levantadas as areas reflo 
restadas do Estado do Mato Grosso do Sul, tendo como base as 
imagens do satelite LANDSAT do ano de 1978 e cujas orbitas/pon 
tos foram 220/25 e 26, 234/25 e 26 e 248/25 e 26. 

Nos arquivos da Delegacia do IBDF em Campo Gran 
de, S, foram levantados os dados contidos nos pro je tos de re 
florestamen^r ao nivel de empresa reflorestadora e propriedade 
rural. Tomando ccmc base crrtografica as cartas da D. S.G. (1973) , 
4 'oram demarcados os limites das fazendas, com a f inalidade de 


/A .COolS 


auxiliar a localizagao e a delimitagao nas iraagens LANDSAT. Fo 
ram feitos sobrevoos e coletadas informagoes no campo para aju 
dar a fase de interpretagao das imagens. 

Nos canais 5 e 7, escalas 1:250.000, os ref lores 
tanentos forammapeados ■’onivel de genero C Pir.ua spp e Eucalyptus 
epp) por propriedade rural. Baseado em informagoes de campo fo 
ram mapeadas as areas preparadas para serem ref lorestadas . 

A quantificagao das classes foi realizada, atra 
ves da tecnica de contagem de pontos, utilizando uma grade mi 
limetrada. 


As areas (ha) do ref lores tamento foram apresenta 
das em modulos correspondentes as imagens que cobrem a area de 
estudo, conforme mostra a Tabela 3. 

A diferenga relativa entre as areas mapeadas nas 
imagens e os dados de verdade terrestre foram da ordem de 2 , 58% . 










5 - LEVANTAMENTO DAS AREAS RE7L0RESTADAS - SEGUNDA FASE 


Nos anos de 1980 e 1981, foram atualizados os da 
dos contidos nos projetos de ref lorestamento na Delegacia Regio 
nal do I.B.D.F. em Campo Grande - MS. Ao nivel de propriedade 
rural, foram mapeadas as areas ref lorestadas utilizando as ima 
gens dos anos de 1979 e 1980, na escala 1:250.000, canais 5 e 
7, e como suporte recorreu-se as imagens dos anos de 1975, 19 76, 
1977 e 1978. Deste modo, aliando-se as informa^oes de campo as 
caracteristicas das areas ref lorestadas das imagens LANDSAT , de 
finiram-se duas classes em funijao do gene.o ( Pinus e Eucalyptus] 
e seis classes em relagao a idade do plantio, conforme mostra 
do na Tabela 4. 

% 

TABELA 4 


CLASSES DE REFLORESTAMENTO EM FUNCAO DO GENERO E IDADE 


CLASSES 

IDADE (ANO) 

OE/OP 

AREA PREPARADA 

1E/1P 

0 a 1 

2E/2P 

1 a 2 

3E/3P 

2 a 3 

4E/4P 

3 a 4 

5E/5P 

4 a 5 

6E/6P 

> 5 



Em relaq:ao a homogeneidade do plantio, estabele 
ceu-se o Indice de Exposi^ao do Terreno (I.E.T.) e def iniram-se 
5 classes conforme mostrado na Tabela 5. 

TABELA 5 

CLASSES DE REFLORESTAMENTO EM FUNC&0 DO I.E.T. 



Indice de expos igAo do terreno 

0 

a 

20% 

20 

a 

40% .. 

40 

a 

60% 

60 

a 

80% 

80 

a 

100% 


As classes de idade de 0 a 1 nao foi associado o 
I.E.T., porque o ref lorestamento ainda nao esta em fase de re 
cobrimento do solo. 

Com a contribui^ao do genero, idade e I.E.T., es 
tabeleceu-se uma legenda com 54 classes de ref lorestamento. 


0 total da area reflorestada em ha, por classe,e 
apresentado na Tabela 6. 
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Area (ha) por classe de idade de reflorestamento 


NO ESTADO DO MATO GROSSO DO SUL 


Area total por 

CLASSE (ha) 


7 . POP CLASSE EM 
re lac ao X Area 
total 



Area total por 
classe de idadf 

(bn) ^ 


14.475, 00 


50.906,25 


59 . 831,25 


7 . DAS CLASSES DE 
IDADE EM FJIUCAO A 

Area total 


44 . 112,50 


73 . 143,75 


66 . 118,75 


46 . 568,75 


6.675.00 


10.262,50 


10 . 968,75 


2 . 575,00 


5 . 931,25 


2 . 337,50 


2 . 806,25 


396 . 912,50 • 



2,76 

4 

0,65 


1,49 


0,64 


0,71 
































































































































6 - multiestAgio 


Atualmente, no Estado do Mato Grosso do Sul,esta 
sendo implementado um procedimento para estimar o volume da ma 
deira denominado Amcstragem em Multiplo com Probabilidade Varia 
vel (Langley, 1975) . 0 objetivo deste metodo e avaliar, ao nl 
vel regional, o volume de madeira, de modo bem rapido. 

A formula geral e o esquema para estimar o volume 
de madeira, utilizando amostragem em tres estagios (dados de sa 
telite, aeronave e campo) , com probabilidade variavel ,conforme 
Langley (1975), sao ilustrados na Figura. 6. 

6.1 - PRIMEIRO ESTAGIO 

No primeiro estagio (nivel orbital), procede-se a 
classificagao das areas ref lorestadas , onde, poster iormente, e 
sobreposta uma grade que permitira uma estimativa de area nas 
Unidades Primarias de Amostragem (UPA) . Supoe-se, por exemplo, 
que se trabalhe com 5 UPAs. Os valores para determinar a proba 
bilidade de selegao sao mostrados na Tabela 7. 


19 ESTAGIO 29 ESTAGIO 39 ESTAGIO 

IMAGEM ORBITAL FOTOGRAFIAS AEREAS CAMPO 

Unidade Primaria de Unidade Secundaria Unidade Terciaria de 

Amostragem (UPA) de Amostragem(USA) Amostragem (UTA) 
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TABELA 7 


VALORES PARA DETERMINAR PROBABILIDADE DE SELECAO DE UPA 



Area (ha) 
( 2 ) 

SOMA 

ACUMULADA 

(3) 

INTERVALO 

(4) 

PROBABILIDADE 
DE SELECAO 
(5) 

36 

36 

1 - 36 

36/135 

25 

61 

37 - 61 

25/135 

15 

76 

62 - 76 

15/135 

43 

119 

77 -119 

43/135 

16 

135 

120 -135 

16/135 


Na Tabela 7, a coluna (11 contem a lista dos nu 
meros em sequencia das UPAs. A coluna (2) mostra os valores de 
areas (ha) de ref lorestamento em cada UPA. A coluna (3) contem 
a f.oma acumulada das areas (ha) da coluna (2). A coluna (4)mos 
tra os valores de intervalos de cada UPA determinados na coluna 
(3) . Na coluna C5) verifica-se que as probabilidades de sele 
9 S 0 , para qualquer UPA, sao variaveis; visto que os seus valo 
res sao determinados pela razao entre a area (ha) , contidos em 
cada UPA, e a area (ha) total das UPAe. Como a porcentagem de 
area e o volume de ref lorestamento, em cada UPA, estao direta 
mente relacionados , verifica-se que, para uma alta porcentagem 
de area reflorestada, existe uma porcentagem alta de volume de 
madeira e, assim, uma alta probabilidade de selecao. 













Os intervalos (coluna 4) sao baseados totalmente • 
nos valores de area (ha) das UPAs. Analisando a UPA numero 1, 
o intervalo de selegao e definido come de 1 a 36. Assim, a pro 
babilidade de selegao para a UPA numero 1 e 36/135, sendo que 
o valor 135 representa a somatSria acumulada na coluna (3) . 

Similarmente, o intervalo para a UPA numero 2 e 
37-61, que representa a probabilidade de selegao de 25/135. 0 

mesmo procedimento e usado para estabelecer os intervalos para 
as outras UPAs. 

Os numeros aleatorios, uniformemente distribui 
dos, sao gerados para cair no intervalo [1, soma acumulada] , on 
de soma acumulada C135 ha) na coluna 3 e o total em (ha) de 
areas ref lorestadas. 

Estes numeros aleatorios sao usados para selecio 
nar as UPAs que sao consideradas no segundo estagio. Com isto, 
e feita a demarcagao dos numeros aleatorios no intervalo [l,so 
ma acumulada] . A UPA selecionada tern a probabilidade de selegao 
proporcional a area (ha) de ref lores tamento e, assim, aprobabi 
lidade de selegao e definida como: 

. _ area (ha) de ref lorestamento para a UPA 
' total da area (ha) de ref lorestamento pa 
ra a area de estudo. 

A area (ha) de ref lorestamento de cada UPA (que 
e diferente para cada UPA) determina a probabilidade de selegao 
e o intervalo para cada UPA. 


0 numtro de UPAs a ser selecionado e definido pe . 
la seguinte formula: 

. t 2 . (CV) 2 

m - , 

(EE) 2 

onde : 

m = numero de UPAs selecionadas no primeiro estagio, 
t = valor "t". Para inventario florestal em extensas areas 
usa-se o valor ”t" = 1 (Wiant Junior, 1976), 

C.V.= coeficiente de variagao, 

•E. E. = erro de amostragem esperado. 

% 

6.2 - SEGUNDO ESTAGIO 


As UPAs selecionadas sao demarcadas em bases car 
tograficas planialtimetricas e subdivididas em quadrados ou re 
tangulos que correspondem as Unidades Secundarias de Amostragem 
(USA). Em fotografias aereas inf ravermelhas coloridas, escala 
1:20.000, as USAs sao delimitadas e. avaliadas em termos de area 
(ha). Os dados assim colecionados e analisados para cada USA, 
dentro de cada UPA, sao selecionados da mesma maneira que no 
primeiro estagio. 


6.3 - TERCEIRO ESTAGIO 


Cada USA selecionada e subdividida em quadradcs 
e retangulos que correspondem as Unidades de Amostragem Tercia 
rias (UTA) . As fotografias aereas inf ravermelhas coloridas, am 


pliadas para a escala de 1:10.000 ou 1:5.000, servem como base 
para a demarcate* das UTAs. A selegao deste estagio e baseada 
no numero ou na densidade de copa das areas ref lorestadas emca 
da UTA. Os dados asslm colecionados sao analisados para cada 
UTA, dentro de cada USA, do mesmo modo que nos estagios ante 
riores. Cada UTA selecionada e demarcada no campo, onde se ob 
tern dados para determinagao de volume de madeira. 

Para estimar o volume total de madeira, utiliza 
-se a formula apresentada na Figura 5. 

Segundo Langley (1975) , para estimar a variancia 
contida na estimativa do volume total (tft) , utili za-se a seguin 
te equagao: 


Var (tft) = — 

m 

E 

Vi 

2 

- m 

m (m-1) 

f — i 

ii 

•H 

1 

_ Pi _ 
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